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Abstract 

In the context of the density functional theory we consider the single particle ex- 
citation spectra of electron systems. As a result, we have related the single particle 
excitations with the eigenvalues of the corresponding Kohn-Sham equations. In a case 
when the exchange correlation functional is approximated by the exchange functional 
the coupling equations are very simple, while the single particle spectra do not coin- 
cide neither with the eigenvalues of the Kohn-Sham equations nor with the ones of the 
Hartree-Fock equations. 

PACS: 31.15.Ew; 31.50.+W Keywords: Density functional theory; Effective interaction; Ex- 
citation spectra 

The density functional theory (DFT), that originated from the pioneer work of Hohen- 
berg and Kohn has been extremely effective in describing the ground state of finite many- 
electron systems. Such a success gave birth to many papers concerned with generalizations 
of DFT, which would permit to describe the excitation spectra also ( see e.g. [2 - 5]). The 
generalizations, on the theoretical grounds, originated mainly from the Runge-Gross the- 
orem, which helped to transform DFT into the time-dependent density functional theory 
TDDFT 0. Both, DFT and TDDFT, are based on the one-to-one correspondence between 
particle densities of the considered systems and external potentials acting upon these par- 
ticles. For the sake of simplicity and definiteness, while dealing with TDDFT, we assume 
that the time dependent part of the considered electronic system's density p(r, t), which is 
created under the action of an external time- dependent field Xvextij^^t), is developing from 
the system's ground state. Unfortunately, the one-to-one correspondence establishes only 
the existence of the functionals in principle, leaving aside a very important question on how 
one can construct them in reality. This is why the successes of DFT and TDDFT strongly 
depend upon the availability of good approximations for the functionals. This shortcoming 
was resolved to large extent in p|, |^, § where an exact equations connecting the action 
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functional, effective interaction and linear response function were derived. But the linear 
response function, containing information of the particle-hole and collective excitations, does 
not directly present information about the single particle spectrum. 

The main goal of this Letter is to derive equations describing single particle excitations 
of multi-electron systems, using as a basis the exact functional equations. As a result, by 
constructing more and more accurate functionals, one will among other things treat more 
and more accurately the single particle excitations. 

Let us briefly outline the equations for exchange-correlation functional Exc[p\ of the 
ground state energy and exchange-correlation functional Aj;c[p] of the action A[p\ in the case 
when the system in question is not perturbed by an external field. In that case an equality 
holds 

E^c[p\ = A^c[p]\pir,u;=0), (1) 

since A^^ is also defined in the static densities domain. The exchange- correlation functional 
E^cIp] is defined by the total energy functional E[p] (see e.g. 0) as 

E[p] ^ To[p] + \l P^p^P^dv.dv, + E,,[p], (2) 

where To[p] is the functional of non- interacting Kohn-Sham particles. The atomic system 
of units e = m = ^= lis used in this paper. The exchange-correlation functional may be 
obtained from: 

\ r 1 dw 

Exc = -7: dg' [x{Yi,T2,iw,g') +2'Kp{vi)5{w)5{vi-T2)]-, -—dTidvi. (3) 

2 JO J |ri —Y2\ 2,71 

Equation (3) presents the well-known expression for the exchange-correlation energy of a 
system ( see e.g. 0, ||, |^), expressed via linear response function x{^i^'^2^iw^g'), with g' 
being the coupling constant. In order to consider (3) as describing E^c, the only thing we 
need, is the ability to calculate the functional derivatives of E^c with respect to the density. 
According to (3), it means an ability to calculate the functional derivatives of the linear 
response function x with respect to the density p(r, tu) which was developed in 0, The 
linear response function is given by the integral equation, 

X(ri,r2,u;,5() = xo(ri,r2,u;) + j xo(ri, r'^, u;)i?(r r 2, t^, fi')x(r2, r2, tu, 5')rfr'i(ir2, (4) 

with xo being the linear response function of non-interacting Kohn-Sham particles, moving 
in the single particle time-independent field [0, §. It is evident that the linear response 
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function x{9 = 1) tends to the linear response function of the system in question as g goes 
to 1. The exact functional equation for R{ri,r2,u!, g) 0, |^ looks as follows, 



R{Ti,V2,uj,g) 



|ri - ral 



(5) 



/ / xir[,r2,iw,g')—^-—dri'dr2'^dg'. 
J Jo \ri — 2tt 



2 5p{ri,uj)6p{r2, -uj) 

Here R{ri,r2,uj, g) is the effective interaction depending on the coupling constant g of the 
Coulomb interaction. The coupling constant g in eq. (5) is varied in the range (0 — 1). The 
single particle potential Vxc, being time-independent, is determined by the following relation 



VrAr 



6p{r[ 



-E 



(6) 



Here the functional derivative is calculated at p = po with po being the equilibrium density. 
By substituting (3) into (6), it can be shown that the single particle potential v^c has the 
proper asymptotic behavior ||, 0, 



1 

r 



(7) 



The potential Vxc determines the energies Si and wave functions 0i of fictitious Kohn-Sham 
particles. 



(r), 



that in turn form the linear response function Xoi 
Xo(ri,r2,u;) = ^ n,(l - nfc)0- (ri)0i(r2)0fc(r2)0yk(ri) 
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1 



u -Uik + irj u + Uik - irj 



and the real density of the system p. 



P(r) = E 



(8) 



, (9) 



(10) 



Here are the occupation numbers, V^xt contains all external single particle potentials of 
the system, say the Coulomb potentials of the nuclei. Then, Vh is the Hartree potential. 



VHir) 



J r — ri 



:iii 
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and ujik is the one-particle excitation energy Uik = Sk 



7] being an infinitesimally small 



positive number. A few remarks regarding the well-known paradox 0, ^ are in order here. 
The linear response function x ci^nd its inverse ^"^^ the noncausal function, with the 
inverse being given by 0, H], 



X ^(ri,r2,cj) 



5^A[p\ 



-uo 



(12) 



These functions can be used on equal footing with the causal ones without any contradictions 
related to the interplay between the causality and symmetry properties of the linear response 
functions 0, |lT], [1^. While the inverse, given by eq. (12), can be used to construct both 
causal and noncausal linear response functions. The described above equations (2-5) solve 
the problem of calculating E^ci the ground state energy and the particle-hole and collective 
excitation spectra of a system without resorting to approximations for E^c based on addi- 
tional and foreign to the considered problem calculations such as Monte Carlo simulations, 
or something of this kind. We note, that using these approximations, one faces difficulties 
in constructing the effective interaction of finite radius and the linear response functions 
0, H . On the base of the suggested approach one can solve these problems. For instance, in 
the case of a homogeneous electron liquid it is possible to determine analytically an efficient 
approximate expression Rrpae for the effective interaction i?, which essentially improves 
the well-known Random Phase Approximation 0, § by taking into account the exchange 
of electrons properly, thus forming the Random Phase Approximation with Exchange. The 
corresponding expression for Rrpae is as follows 



R 



where 



RPAE 



gn 



RE{q,9,p) = 2 



Pf 



In 



^ + RE{q,g,p), 



12pi 



4p 



2pF_ 

3g 



In 



2pF - q 



2pF + q 



1 

+ 3 



(13) 



(14) 



Here the electron density p is connected to the Fermi momentum by the ordinary relation 
p = pp/STr"^. Having at hand the effective interaction RppAEiq, 9, p), one can calculate the 
correlation energy Ec per electron of the electron gas with the density Vg. The dimensionless 
parameter = ro/ap is usually introduced to characterize the density, with tq being the 
average distance between electrons, and is the Bohr radius. The density is high, when 
r, < 1. 
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Table 1. Correlation energy per electron in eV of an electron gas of the density r^. The 

are compared with the RPA calculations and with 



Monte Carlo results [13 



the results of Srpa denotes the results of RPA calculations, and Srpae 
denotes the results of the calculations when the effective interaction R was 
approximated by Rrpae 0- 







£rpa 


^RPAE 


1 


-1.62 


-2.14 


-1.62 


3 


-1.01 


-1.44 


-1.02 


5 


-0.77 


-1.16 


-0.80 


10 


-0.51 


-0.84 


-0.56 


20 


-0.31 


-0.58 


-0.38 


50 


-0.16 


-0.35 


-0.22 



As can be seen from Table 1, the effective interaction Rrpae{<1i p) permits to describe the 
electron gas correlation energy in an extremely broad interval of density variation. Note, 
that even at = 10 the mistake is no more than 10% as compared to Monte Carlo calcula- 
tions, while the result becomes almost exact at = 1 and is absolutely exact when — > 

i- 

Now let us calculate the single particle energies e^, that, generally speaking, do not 
coincide with the eigenvalues Si of eq. (8). Note that these eigenvalues do not have the 
physical meaning of a particle energy (see e.g. []T|). To calculate the single particle energies 
one can use the Landau equation [0 

5E , , 

J- = e^. (15) 
drii 

Eq. (15) can be used since, as it follows from eqs. (4,9,10), the linear response functions x 
and Xo; and the density p depend upon the occupation numbers. Thus, one can consider the 



ground state energy as a functional of the density and the occupation numbers 



E[p{r),ni] = Tfc[p(r),n,] + VH{r)p{r)d{r) + J Kxt(r)p(r)d(r) + E,,[p(r), n^]. (16) 

Here is the functional of the kinetic energy of non-interacting Kohn-Sham particles. Upon 
substituting eq. (16) into (15) and using eqs. (3,8), one gets, that the single particle energies 
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ei can be presented by the following expression, 



< (Pi\VxcWi > -7:^ 

2 orij 



X(ri, T2, iw) + 27rp(ri)^(w)^(ri - 



dwdg'dridr2 



27r 

(17) 

The first and second terms on the right in eq. (17) is linked with the derivative of func- 
tional Tfc with respect to occupation numbers n^. To calculate the derivative we consider an 
auxiliary system of non-interacting particles moving in a field U{r). Then its ground state 
energy Eq is given by the following equation 

Eo = Tk + J U{r)p{r)dr. (18) 

Varying Eq with respect to the occupation numbers one gets the desirable result, 

= ^ = ^+ < m\<l>^ >, (19) 
OTli OUi 

provided U = Vh + V^xt + Vxc- The third term is related to the contribution coming from 
Er^c defined by eq. (3). The derivative Sx/Srii, as it follows from eq. (4), is given by, 

^ = ^ + -^Rx + Xo^X + XoR^- (20) 
orii orii OUi orii orii 

In eq. (20) for the sake of brevity we omit the spatial integrations. The variational derivative 
Sxo/Srii is calculated directly from (9) and has the simple functional form, 

'^^"^'^''"""'^ = [G'o(ri, T2, 00 + 6,) + Go(ri, r^, -u + e,)] 0*(ri)0.(r2), (21) 

with Go(ri,r2,a;) being the Green function of non- interacting electrons moving in the 
single particle potential Vh + v^c + Kxt- We do believe, and some evidences can be found 
in [0, that the contribution coming from 6R/6ni is small. Nonetheless even having omitted 
this term, we have still to deal with rather complex eq. (17). Therefore, it is of interest to 
illustrate the general consideration of the one-electron energies with a simple and important 
example, when only the exchange part E^ of the total exchange-correlation functional is 
selected to be treated rigorously, using an approximation for Ec = E^c — E^. Taking into 
account calculations of given in Table 1, and the local density approximation, one gets, 

EM = J P(r)ee((r))rfr. (22) 
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For the functional we have an exact expression [|], |T6 

If 1 dw 

E4p] = -7; / [Xo(ri,r2,iw) + 27rp(ri)5(w)5(ri - ra)] ^ -—dridr2. 

2 J |ri — r2| 

As it is seen from eq. (6), single particle potential v^c takes the form, 

f^c(r) = Wx(r) +fc(r), 
with the potentials Vx{r) and fc(r) given by the following expressions 

6E^ _ 5Er. 



6p{r) 



6piv) 



(23) 



(24) 



(25) 



The potential fx(i") can be calculated exactly ||, |T6|, |T3, while there are quite suitable 
approximations to the potential fc(r) (see e.g. P). Presenting single particle potential 
fa;c(r) in such a way as (25), we simplify the calculations a lot, preserving at the same time 
the asymptotic condition (7), namely Vxdf' — ^ 00) — 1/r. This condition is of crucial 
importance for calculations of Ei related to the Kohn-Sham vacant states (see eq. (8)). In 
the same way as eq. (17) was derived, one gets. 



Ei- < (j)i\v^\ 



> - 



1_5_ 

2 6ni 



Xo(ri, r2, iw) + 27rp{ri)6{w)6{ri - r2) 



ri - r2 



dw dridr2 
2^ ' 



(26) 



After straightforward calculations one arrives at rather simple result for the single particle 
spectra, that are to be compared with experimental results. 



< 4>i\Vx\ 



> 



[ri)0^(r2)0^(r2)0fc(ri; 
ki - r2| 



dridr2. 



(27) 



single particle levels e^, given by eq. (27), resemble the ones that are obtained within the 
Hartree-Fock (HF) approximation. If the wave functions (pi would be solutions of the HF 
equations and the correlation potential fc(r) would be omitted, the energies would exactly 
coincide with their eigenvalues Ei. But this is not the case, since 0, are solutions of eq. (8), 
and the energies do not coincide with either HF eigenvalues or with the ones of eq. (8). 

In summary, we have shown that it is possible to calculate the single particle excitations 
within the framework of DFT. The developed equations allow for calculations of the single 
particle excitation spectra of any multielectron system such as atoms, molecules and clusters. 
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We anticipate as well that these equations when applied to solids will produce a quite reason- 
able results for the single particle spectra and energy gap at various high-symmetry points 
in the Brillouin zone. We have also related the eigenvalues of the single particle Kohn-Sham 
equations to the real single particle spectrum. In the most straightforward case, when the 
exchange functional is treated rigorously while the correlation functional is taken in the local 
density approximation, the coupling equations arc very simple. The single particle spectra 
do not coincide either with the eigenvalues of the Kohn-Sham equations or with the ones 
of the Hartree-Fock equations, even provided the contribution coming from the correlation 
functional is omitted. 

This research was funded in part by INTAS under Grant No. INTAS-OPEN 97-603. 
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